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1.  INTRODUCTION 


The  final  product  of  powder  metallurgy  or  ceramics  processing  techniques  is  typically  a  solid  made 
up  of  powder  particles  consolidated  by  hot-pressing  or  sintering.  The  mechanical  properties  of  the  product 
are  often  highly  dependent  on  the  quality  of  intergrain  bonding  which,  in  turn,  may  depend  strongly  on 
near-surface  impurities  present  on  the  powder  particles.  Thus,  the  presence  of  a  thin  oxide  or  metal 
contaminant  layer  may  have  a  drastic  effect  on  the  outcome  of  a  powder  consolidation  process.  There 
exist  techniques,  based  on  the  scattering  of  high-energy  ions,  that  are  capable  of  characterizing  such 
contaminant  layers  on  flat  surfaces.  The  ultimate  goal  of  this  work  is  to  adapt  these  methods  so  that  they 
can  be  used  on  powder  particle  surfaces. 

When  a  beam  of  high-energy  ions  strikes  the  surface  of  a  solid  sample,  most  of  the  ions  simply  lose 
energy  through  numerous  small  interactions  with  target  electrons  and  eventually  come  to  a  halt  some 
distance  below  the  surface.  A  critical  few,  however,  will  scatter  off  a  target  nucleus  and  out  of  the  sample 
where  they  can  be  detected.  The  shape  of  the  energy  spectrum  of  these  scattered  particles  can  often  be 
interpreted  so  as  to  yield  information  concerning  the  distribution  of  elements  near  the  surface  of  the 
sample.  In  many  such  "ion  beam  analysis"  techniques,  the  interpretation  is  carried  out  through  computer 
simulation  of  the  particle-solid  interactions  that  take  place.  Depending  on  the  type  of  technique  used,  the 
intemuclear  scattering  events  may  be  due  to  simple  electrostatic  repulsion  (as  is  the  case  for  Rutherford 
Backscattering  Spectroscopy  [RBS]  [Chu,  Mayer,  and  Nicolet  1978])  or  due  to  nuclear  reactions  (Nuclear 
Reaction  Spectroscopy  [NRS]  [Niiler  and  Kecskds  1989]).  These  techniques  utilize  ions  with  energies  in 
the  megaelectronvolt  (MeV)  range,  resulting  in  useful  detection  ranges  on  the  order  of  a  few  microns  from 
the  surface. 

While  ion  beam  techniques  such  as  RBS  and  NRS  can  yield  considerable  information  on  a  sample’s 
surface  condition,  current  methods  of  ion  beam  analysis  generally  require  a  flat,  laterally  uniform  surface. 
This  is  because  the  simulation  process  becomes  extremely  complicated  when  a  surface  with  something 
other  than  this  simplest  of  configurations  is  considered.  This  report  addresses  methods  of,  and  issues 
involved  with,  adapting  ion  beam  analysis  techniques  for  use  on  spherical  and  cylindrical  particle  surfaces. 

Two  systems  for  dealing  with  curved  surfaces  will  be  discussed;  both  attempt  to  reduce  the  problem 
of  simulating  backscattered  particle  energy  spectra  arising  from  spheres  and  cylinders  to  a  point  where  flat 
surface  simulations  can  be  used  to  construct  them.  The  objectives  of  both  are  to  extract  information  on 
the  near-surface  radial  concentration  profiles  of  these  surfaces  by  fitting  experimentally  obtained  spectra. 
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The  first  method  (which  is  similar  to  one  used  previously  [Niiler  and  Kecsk6s  1989])  straightforwardly 
approximates  the  curved  surface  with  a  finite  array  of  flat  surfaces,  each  tilted  in  a  particular  direction 
relative  to  the  beam  direction.  The  distribution  of  "tilts"  is  characteristic  of  the  type  of  surface  in  question 
(i.e.,  sphere  or  cylinder).  A  flat  surface  simulation  is  generated  for  each  of  these  surfaces  and  these  are 
summed  to  obtain  the  final  result.  This  method  will  be  referred  to  as  the  "polyhedron  approximation." 
The  procedure  described  is  a  typical  "trial  and  error"  method  where  a  radial  profile  is  first  chosen  and  a 
spectrum  representing  it  is  then  simulated  and  compared  to  actual  data.  If  the  fit  is  imperfect,  then  a  new 
radial  profile  is  chosen,  another  simulation  is  performed,  and  so  on,  until  an  acceptable  fit  is  achieved. 

The  second  set  of  methods  arose  from  the  observation  that  a  spectrum  obtained  from  a  curved  surface 
could  often  be  successfully  fit  under  the  (incorrect)  assumption  that  the  surface  was  flat  (Chu,  Mayer,  and 
Nicolet  1978).  This  suggests  the  existence  of  an  "equivalent"  flat  surface  depth  profile  that  mimics  the 
results  for  a  given  radial  profile  on  a  sphere  or  cylinder.  The  "flat  model"  encompasses  two  types  of 
procedures  that  make  use  of  this  idea.  The  first  procedure  is  based  on  the  construction  of  an  equivalent 
or  nearly  equivalent  flat  surface  depth  profile  from  an  arbitrary  radial  profile.  Only  a  single  flat  surface 
simulation  need  then  be  performed  in  order  to  obtain  the  final  spectrum  representative  of  that  radial 
profile.  As  in  the  polyhedron  approximation  procedure,  a  spectrum  simulated  this  way  can  then  be 
compared  to  actual  data  so  as  to  judge  the  suitability  of  the  radial  profile  chosen.  The  second  "flat  model" 
procedure,  which  is  only  applicable  in  certain  cases,  comes  at  the  problem  from  the  other  direction.  In 
it,  a  spectrum  from  a  curved  surface  is  immediately  fit  as  if  it  were  actually  from  a  flat  surface.  The 
desired  radial  profile  functions  are  then  derived  from  the  "false"  flat  surface  profile  obtained  in  the  fit, 
using  a  simple  mathematical  procedure.  In  this  situation,  no  extraordinary  simulation  tools  are  needed  at 
all.  These  "flat  model"  methods  will  be  compared  to  the  polyhedron  approximation  method,  and  the 
limitations  of  each  will  be  discussed. 

While  the  methods  discussed  in  this  report  can  be  applied  to  a  variety  of  ion-beam  techniques, 
examples  shown  here  will  be  confined  to  2-MeV  4He  RBS.  Spectra  simulations  will  be  constructed  with 
the  aid  of  the  flat  surface  RBS  simulation  program  RUMP  (Doolittle  1985).  It  will  always  be  assumed 
that  the  spheres  or  cylinders  to  be  studied  are  completely  spanned  by  the  probing  ion  beam;  the  use  of 
microbeams  is  not  considered.  It  will  also  be  assumed  that  surface  layers  are  either  spherically  or 
cylindrically  symmetric.  While  the  assumptions  of  large  radii  samples  and/or  180°  backscattering  will 
often  be  invoked,  either  because  they  are  required  or  convenient,  alternate  situations  will  also  be  discussed. 
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2.  POLYHEDRON  APPROXIMATION 


Currently,  backscattering  spectrum  simulations  are  generally  restricted  to  the  simplest  sample 
geometry:  a  flat,  laterally  uniform  surface.  Typically,  the  only  geometrical  modification  allowed  is  a 
change  in  the  beam  to  surface  normal  angle,  or  tilt.  In  the  case  of  powder  surface  analysis,  the  sample 
surface  is  better  described  by  a  layer  of  close-packed  spheres  where  each  sphere  presents  a  range  of 
surface  normal  angles  to  the  beam.  In  such  cases,  the  effective  thickness  of  a  uniform  layer  will  vary 
across  the  spherical  surface.  Consequently,  if  a  thin  layer  on  a  flat  surface  results  in  a  peaklike  feature 
in  a  backscattering  spectrum,  the  same  layer  on  a  sphere  will  result  in  a  broader  peak  with  a  long,  low- 
energy  tail  (Niiler  and  Kecskds  1989).  The  method  to  be  described  here  accounts  for  this  by 
approximating  a  sphere  (or  cylinder)  with  an  array  of  separate  flat  surfaces,  each  having  an  appropriate 
beam  to  surface  normal  angle  and  size.  This  is  somewhat  similar  to  the  methods  mentioned  in  Niiler  and 
Kecsk6s  (1989)  and  Feldman  and  Zinke-Allmang  (1990).  Since  the  distribution  of  surface  normals  is  the 
same  for  all  spheres,  only  one  sphere  need  be  considered.  In  the  case  of  aggregates,  it  will  be  assumed 
that  the  ion  beam  strikes  the  sample  obliquely,  so  the  shadowing  of  one  sphere  by  another  is  not 
considered. 

The  following  describes  one  way  of  arriving  at  an  appropriate  array  of  flat  surfaces.  This  is 
accomplished  by  dividing  the  surface  of  a  sphere  into  small  regions  and  approximating  each  region  with 
a  flat  surface  having  the  same  mean  beam-to-normal  and  detector-to-normal  angles  as  well  as  the  same 
frontal  area  exposed  to  the  beam.  Due  to  the  symmetry  of  the  situation,  only  one  quarter  of  a  spherical 
surface  need  be  described:  we  choose  the  region  0<z<R,  0<y<R,  — R  <  x  <  R  ,  where  R  is  the 
radius  of  the  sphere,  as  illustrated  in  Figure  1.  If  the  beam  travels  in  the  -z  direction  and  the  detector  lies 
in  the  x-z  plane  at  an  angle  of  8  with  the  z-axis,  then  the  beam  to  normal  angle  (j)bn  and  the  detector  to 
normal  angle  <j)dn  at  an  arbitrary  point  on  the  sphere  (x,  y,  z)  are  given  by: 


cos  (j>dn  =  JLsin  8  +  -Lcos8.  (2) 

■ran  R  R 
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One  could  partition  the  surface  of  the  sphere  into  regions  described  by  circles  separated  by  constant 
polar  and  azimuthal  angles  (i.e.,  equally  spaced  "longitude"  and  "latitude"  lines).  However,  this  proves 
to  be  relatively  inefficient  for  the  proposes  of  simulating  backscattering  spectra.  This  is  due  to  the  fact 
that  the  effective  thickness  of  a  uniform  layer  varies  slowly  near  the  center  of  the  sphere  and  quickly  near 
the  edge  of  the  sphere.  In  other  words,  much  of  the  structure  seen  in  a  simulated  spectrum  arises  from 
regions  near  the  edge  of  the  sphere  and  thus  it  is  more  efficient  to  have  a  high  density  of  regions  near  the 
edge  and  a  low  density  near  the  center.  One  way  of  accomplishing  this  is  to  choose  regions  described 
by  circles  separated  by  constant  z,  rather  than  circles  separated  by  constant  polar  angle.  We  choose  to 
divide  the  quarter  sphere  into  2n2  regions,  with  n  divisions  in  the  z  direction  and  2n  divisions  in  azimuthal 
angle.  The  midpoints  of  these  regions  have  z  and  x  coordinates  given  by: 


(.  1  Yr  "I 

Z;  =1  -  _  - 

1  I  2Anj 


,  i  =  l,...,n 


(3) 


=  /R2  - 


COS 


2 


iV*'' 


v2n  / 


,  j  —  l,...,2n. 


(4) 


Here  the  index  "i"  refers  to  the  division  in  z  and  "j”  refers  to  the  division  in  azimuthal  angle.  Under  the 
assumption  of  a  uniform  beam  completely  spanning  the  sphere,  the  total  charge  deposited  in  each  region 
will  be  proportional  to  the  its  frontal  area  exposed  to  the  beam.  Here  the  fraction  of  the  total  area  exposed 
by  each  region  is  given  by: 


=  (2i  -  1) 


fi  1 


(5) 


This  quantity  does  not  depend  on  the  index  j.  Thus  the  final  spectrum  is  simulated  by  summing  2n  flat 
surface  simulations  each  with  angles  given  by  Equations  1-4  and  with  the  fraction  of  the  total  charge 
deposited  given  by  Equation  5.  Note  that  the  dependence  on  R  drops  out.  This  approximation  woiks  best 
when  a  large  number  of  segments  are  used;  we  find  that  50  to  100  separate  simulations  will  result  in  a 
level  of  detail  sufficient  for  fitting  purposes.  In  the  particular  case  of  180°  scattering  (8  =  0),  azimuthal 
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symmetry  makes  much  of  the  process  redundant,  so  that  only  n  calculations  need  be  performed,  with 
parameters: 


V 


1) 

V 

2 

n 

aj  =  (2i  -  1) 


r.i 


(7) 


A  similar  process  can  be  built  up  for  cylindrical  surfaces — if  the  axis  of  the  cylinder  corresponds  to 
the  y-axis,  then  the  salient  differences  are  that: 


1  - 


z 


,  j  — 1,1 


(8) 


ai  = 


1  - 


i  -  1 


\2 


\  n  J 


1  - 


(\  'V 


,  1  l,..;,n. 


(9) 


While  the  above  method  may  resemble  the  approximation  of  a  sphere  or  a  cylinder  with  an  irregular 
polyhedron,  it  differs  in  one  major  respect:  it  in  no  way  accounts  for  the  effects  of  incident  particles 
entering  through  one  facet  and  the  resulting  scattered  particles  leaving  through  another  facet  The  implicit 
assumption  that  these  effects  are  negligible  is  tantamount  to  the  assumption  that  the  sphere’s  radius  is 
much  larger  than  the  probing  range  of  the  ion  beam  technique.  This  restriction  to  large  radii  further 
implies  that  regions  out  of  direct  view  from  the  detector  (<J>dn  >  rc/2)  should  be  ignored.  While  the  total 
amount  of  frontal  area  represented  by  these  inaccessible  regions  may  be  small  (<1%  for  5  =  10°),  they 
may  represent  a  noticeable  portion  of  the  edge  regions  responsible  for  the  low-energy  tails  seen  in 
backscattering  spectra.  Thus,  increasing  the  detector  angle  8  tends  to  decrease  the  yield  in  these  tails. 
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The  RBS  simulation  program  RUMP  (Doolittle  1985)  is  easily  adapted  for  use  with  this  approximation 
because  its  facility  for  stringing  together  subcommands  into  a  macrocommand  obviates  the  modification 
of  its  source  code.  One  need  only  set  up  the  desired  radial  profile  description  in  RUMP’s  simulation 
subprocessor  (as  if  it  were  a  flat  surface  profile)  and  then  invoke  a  macrocommand  that  performs  all  of 
the  separate  simulations  and  sums  them  to  obtain  the  result  for  a  curved  surface.  While  this 
macrocommand  may  be  long  and  complex,  it  can  quickly  be  set  up  in  a  spreadsheet  program.  Figures  2a 
and  2b  demonstrates  the  use  of  this  system  in  conjunction  with  2-MeV  4He+  RBS:  they  show 
experimental  and  simulated  spectra  from  a  0.1 -in-diameter  copper  wire  onto  which  a  Ta  layer  had  been 
ion-plated,  along  with  the  derived  radial  elemental  profiles.  Most  of  the  low-energy  tail  seen  on  the  Ta 
peak  is  due  to  the  cylindrical  shape  of  the  substrate,  although  some  intermixing  is  evident  It  can  also  be 
seen  that  the  Ta  layer  is  heavily  contaminated  by  C,  N,  and  0;  the  proportions  of  which  were  determined 
using  1-MeV  d+  NRS  and  the  simulation  program  PROFILE  (Niiler  and  Kecsk6s  1989). 

The  polyhedron  approximation  has  proven  to  be  a  versatile  method  for  adapting  ion  beam  techniques 
to  the  study  of  contamination  of  curved  surfaces.  It  has,  however,  one  fundamental  restriction  and  one 
practical  limitation:  it  is  restricted  to  use  on  samples  with  large  radii  and  it  tends  to  be  very  slow.  The 
simulation  process  is  slow  because  it  is  often  necessary  to  use  large  numbers  of  separate  flat  surface 
simulations  to  describe  all  of  the  features  seen  in  spectra  adequately,  particularly  the  low-energy  tails  seen 
on  peaks.  The  simulation  seen  in  Figure  2  is  the  sum  of  80  separate  flat  surface  simulations,  and,  while 
RUMP  is  very  fast,  it  required  over  20  min  on  an  IBM-AT  compatible  computer  with  a  12-MHz  80286 
microprocessor  and  80287  math  coprocessor  to  perform  all  of  these.  These  limitations  were  the  motivation 
behind  development  of  the  next  method  to  be  discussed. 

3.  FLAT  MODEL 

3.1  Introduction.  As  seen  in  Figure  2a,  spectra  arising  from  samples  with  rounded  surfaces  tend  to 
resemble  those  from  flat  surfaces  in  which  there  is  intermixing  between  the  surface  and  bulk  regions  (Chu 
and  Mayer  1978).  In  fact,  such  spectra  can  often  be  fit  if  a  flat  surface  is  assumed,  suggesting  that  a 
particular  flat  surface  depth  profile  exists  that  mimics  the  round  surface  case.  This  is  akin  to  the  case  of 
a  flat  surface  with  an  impurity  depth  profile  given  by  a  function  N(f)  and  surface  normal  tilted  at  an  angle 
6  to  the  beam:  if  the  detector  is  at  180°,  then  an  untilted  flat  sample  with  an  impurity  profile  given  by 
N(Ccos0)  will  result  in  the  same  backscattering  spectrum.  The  purpose  of  this  model  is  to  construct  just 
such  an  "equivalent"  flat  surface  profile  that  will  mimic  a  spectrum  from  a  sphere  or  cylinder. 
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The  discussion  that  follows  closely  resembles  that  found  in  Chu,  Mayer,  and  Nicolet  (1978),  except 
that  here  the  sample  may  have  many  constituents  and  the  energy  loss  per  unit  length  (dE/dx)  needn’t  be 
independent  of  energy.  We  begin  by  assuming  180°  backscattering,  so  that  an  incoming  particle  and  the 
resultant  scattered  particle  must  travel  the  same  path  in  the  sample.  Next,  for  the  sake  of  argument,  we 
assume  that  the  surface  and  bulk  values  of  dE/dx  are  the  same  at  any  given  energy  (i.e.,  dE/dx  is 
independent  of  variations  in  composition).  If  this  is  the  case,  then  particles  scattered  from  nuclei  of  type 
"i"  and  having  the  same  energy  "E"  upon  exiting,  the  sample  must  have  scattered  from  the  same  depth 
T  below  the  spherical  surface.  Thus  it  can  be  seen  that  the  backscattering  yield  from  element  "i" 
between  energies  E  and  E  +  dE  will  be  proportional  to  the  number  of  atoms  of  type  ”i"  found  between 
the  corresponding  surfaces  defined  by  the  depths  fi  and  fi  +  dfi.  A  key  point  is  that  die  actual  position 
of  each  of  these  atoms  need  not  be  known,  other  than  that  they  lie  between  the  depths  fi  and  fi  +  dfi. 
Thus,  only  the  average  concentrations  as  a  function  of  depth  are  crucial,  and  this  fact  suggests  that  these 
average  composition  functions  might  serve  as  "equivalent"  flat  surface  profiles  for  the  purpose  of  spectra 
simulation.  Note  that  determining  the  average  composition  vs.  depth  in  a  sphere  given  a  particular  set  of 
radial  elemental  distributions  requires  only  a  geometrical  analysis.  This  does  not  reference  a  specific  type 
of  ion  beam  analysis,  as  in  Chu,  Mayer,  and  Nicolet  (1978),  Shorin  and  Sosnin  (1991),  Takacs,  Ditroi, 
and  Mahunka  (1989),  and  Jex,  Hill,  and  Mangleson  (1990),  although  we  will  confine  examples  to  RBS 
here. 

It  should  be  stressed  that  we  do  not  expect  that  this  method  will  deliver  exactiy  equivalent  flat  surface 
profiles  in  all  cases,  but  rather  it  is  expected  to  yield  useful,  nearly  equivalent  flat  surface  profiles  in  a 
variety  of  cases.  Comparisons  with  the  polyhedron  approximation  will  indicate  when  the  "flat  model" 
method  is  inappropriate. 

We  have  found  that  average  concentration  functions  can  in  fact  be  used  as  equivalent  or  near 
equivalent  flat  surface  depth  profiles,  within  certain  limitations.  The  practicability  of  this  method  is 
limited  by  the  two  assumptions  listed  above.  The  first  requirement,  that  of  180°  scattering,  can  nearly  be 
met  through  the  use  of  an  annular  detector.  The  second  requirement,  that  the  value  of  dE/dx  at  any  given 
energy  be  exactly  the  same  for  all  compositions,  cannot  generally  be  met,  however.  The  question  then 
becomes:  How  similar  must  the  surface  and  bulk  values  of  dE/dx  be  in  order  to  maintain  the  practical 
applicability  of  this  method?  We  have  probed  the  extent  of  this  requirement  by  simulating  spectra  for 
uniform  layers  of  various  elements  on  spheres  of  other  elements  and  then  comparing  these  to  simulations 
generated  using  the  polyhedron  approximation.  As  the  shapes  of  dE/dx  vs.  energy  curves  tend  to  be 
similar,  we  have  arbitrarily  chosen  to  quantify  differences  between  substances  by  comparing  values  of 
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dE/dx  at  2  MeV  only.  Energy  loss  values  are  taken  from  Ziegler,  Biersack,  and  Littmaik  (1985). 
Remarkably,  we  have  found  that  (for  4He  RBS)  the  surface  and  bulk  values  of  dE/dx  must  differ  by  more 
than  a  factor  of  2  before  the  performance  of  this  method  is  appreciably  degraded.  This  allows  for  a  wide 
variety  of  sample  configurations. 

Two  different  types  of  procedures  that  utilize  the  "equivalent"  flat  surface  profile  concept  will  be 
described.  The  first  is  similar  to  the  polyhedron  approximation  in  that  a  radial  elemental  profile  is  first 
assumed  and  a  simulated  spectrum  is  then  generated  for  comparison  with  an  actual  spectrum.  The 
simulated  spectrum  is  generated  by  constructing  a  flat  surface  profile  that  is  "equivalent"  to  the  round 
surface  radial  profile  chosen  and  then  inserting  this  into  a  flat  surface  simulation  program  (e.g.,  RUMP). 
The  second  procedure,  which  is  applicable  in  only  a  relatively  narrow  range  of  situations,  is  even  simpler 
than  the  first  Here  the  spectrum  from  a  round  surface  is  immediately  fit  under  the  (incorrect)  assumption 
that  the  surface  is  flat,  using  standard  flat  surface  simulation  techniques.  The  desired  radial  elemental 
profiles  can  then  be  derived  directly  from  the  "fitted"  (as  opposed  to  "constructed")  equivalent  flat  surface 
profile.  Thus,  no  prior  assumptions  concerning  the  radial  profiles  need  be  made  and  no  extra  steps  are 
inserted  into  the  flat  surface  simulation  process. 

3.2  Constructing  an  "Equivalent"  Flat  Surface  Profile:  N(t).  The  following  describes  a  method  for 
finding  an  "equivalent"  flat  surface  profile  function  N({)  for  the  i*  constituent  of  a  sample,  given  that  its 
radial  profile  is  described  by  the  function  Nj(r).  As  pointed  out  in  the  previous  section,  this  amounts  to 
a  simple  spatial  averaging  procedure. 

As  illustrated  in  Figure  3,  the  focus  of  points  situated  at  a  constant  depth  T  in  a  sphere  (as  measured 
in  the  direction  of  an  incident  ion  beam)  is  a  section  of  spherical  surface  with  the  same  radius  R,  but 
centered  a  distance  (  below  the  center  of  the  original.  We  wish  to  determine  the  average  atomic 
concentration  of  the  i*  constituent  in  the  volume  contained  between  the  surfaces  at  l  and  t  +  dl,  given 
that  its  radial  atomic  density  profile  is  described  by  the  function  Nj(r),  where  r  is  the  radial  distance  as 
measured  from  the  surface.  We  will  restrict  the  calculation  to  spherically  symmetric  elemental  profiles. 
The  total  number  of  atoms  of  type  "i"  found  in  a  narrow  ring-shaped  section  of  this  shell  located  between 
polar  angles  of  0'  and  0'  +  d0'  subtended  at  the  center  of  the  shell  is: 


dNfng  =  Nj(r)27cR2sin  0'  cos  0'  d0'  dfi  , 


(10) 
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Utilizing  Equation  11,  this  can  be  reduced  so  that  the  integral  is  in  terms  of  n 


Ni(i,  R)  =  Jo4Ni  (r)  Fs  (UR)  dr, 


where  the  function  Fs  (r,  t,  R)  is  given  by: 


Fs  (r,  t  R)  = 


-  3Rr2  +  (2R  2 


-  «2)  r  +  R{2 


). 


(14) 


Equation  13  embodies  the  most  crucial  step  of  this  flat  model  procedure.  It  shows  that,  in  order  to 
construct  an  "equivalent"  flat  surface  profile  Nj({),  one  need  only  integrate  the  multiple  of  the  assumed 
radial  profile  function  Nj(r)  and  a  third-degree  polynomial  in  r  (Equation  14).  In  the  extreme  case  where 
the  radius  of  the  sphere  is  much  larger  than  the  probing  range  (i.e.,  R  »  t),  the  function  Fs  reduces  to 
an  even  simpler  form: 


lim  Fs(r,  «,  R)  =  -1  r.  (15) 

R-»~  l2 


Note  that  this  does  not  revert  back  to  a  case  of  a  flat  surface  as  R  — >  because  of  the  implicit 

assumption  that  the  entire  sphere  is  illuminated  by  the  ion  beam  at  all  times.  Thus,  the  entire  range  of 
"apparent"  thicknesses  of  a  uniform  layer  on  the  sphere  is  always  represented. 

A  similar  process  can  be  developed  for  cylinders,  resulting  in  an  integral  transform  as  found  in 
Equation  13,  except  that  the  function  F  for  a  cylinder  has  die  unwieldy  form: 
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When  R  » 1,  this  reduces  to  the  more  convenient: 


lim 

R  ->oo 


Fc  (r,  I,  R)  = 


(r/t2) 

{l  -  (r/4)2 


(17) 


In  order  to  demonstrate  this  procedure,  consider  a  uniform  layer  of  a  material  "1"  on  a  sphere 
consisting  of  a  material  "2,"  so  that: 


Nj(r)  = 


[0 


0  <  r  <  t 
r  >  t 


N2(r)  = 


0  0  <  r  <  t 

N2  r  >  t 


(18) 


where  t  is  the  thickness  of  the  layer  and  Nj  and  N2  are  the  atomic  densities  of  the  two  materials.  The 
"equivalent"  flat  profile  for  the  surface  material  will  then  be  given  by  Equation  13,  the  result  of  which 
is: 
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In  the  R  »  t  limit,  this  reduces  to 


Njfl)  = 


Ni 


N, 


0  <  «  <  t 


(20) 


t  >  t 


Figure  4  illustrates  these  and  other  results  of  this  process  by  comparing  the  shapes  of  several  radial 
profile  functions  N(r)  with  the  shapes  of  their  resultant  "equivalent"  flat  surface  profile  functions  Nj(f). 
All  assume  that  R  » 1.  Examples  of  radial  profile  functions  shown  include  a  uniform  layer  of  thickness 
t  on  a  sphere  and  cylinder,  a  linear  gradient  on  a  sphere,  and  a  complementary  error  function  profile  on 
a  sphere.  The  first  two  cases  resulted  in  analytical  solutions,  while  the  error  function  case  was  solved 
numerically. 

Unlike  the  polyhedron  approximation,  this  method  is  not  restricted  to  large  spheres.  Figure  5  shows 
the  average  density  as  a  function  of  depth  for  a  uniform  layer  on  a  sphere  for  various  ratios  of  radius  R 
to  film  thickness  t.  It  is  apparent  that  the  smaller  the  radius  of  the  sphere,  the  more  prominent  the  tail  of 
the  "equivalent"  profile.  Figure  5  shows  an  extreme  case  (R/t  =  5)  where  the  entire  sphere  would  be 
profiled.  It  is  unlikely,  however,  that  this  model  would  be  useful  for  studying  aggregates  of  particles  with 
diameters  on  the  order  of  the  probing  range,  as  it  does  not  account  for  beam  particles  interacting  with 
more  than  one  sphere  at  a  time.  For  samples  that  contain  particles  with  a  known  distribution  of  radii 
(between  the  limits  R  =  6  and  R  »  i),  it  may  simply  be  necessary  to  average  the  "equivalent"  flat 
surface  profiles  in  order  to  obtain  a  single  "equivalent"  flat  profile  representative  of  the  entire  sample. 
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The  next  step  is  to  actually  simulate  spectra  so  that  they  may  be  compared  to  simulations  based  on 
the  polyhedron  approximation.  In  order  to  insert  the  above  information  into  a  spectrum  simulation 
program  such  as  RUMP  or  PROFILE,  it  is  first  necessary  to  determine  the  fractional  concentration  of  each 
constituent  as  a  function  of  depth: 


fid)  = 


Njd) 

ENjd) 


(21) 


as  well  as  the  total  areal  density  as  a  function  of  depth: 

Ad)  -  J^ENjCOdl.  (22) 

i 

For  a  meaningful  comparison  with  the  polyhedron  approximation  method,  it  is  necessary  to  restrict  all 
cases  to  R  »  {.  For  the  example  of  a  uniform  layer  on  a  sphere  as  described  in  Equations  18  and  20, 
this  leads  to: 


and 


fid)  = 


0  <  I  <  t 
t  >  t 


1  + 
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N, 


_  -  1 
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Njl  0  <  (  <  t 
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The  above  expressions  and  the  program  RUMP  were  used  to  create  spectra  simulations  that  could  be 
compared  to  simulations  generated  using  the  polyhedron  approximation.  A  variety  of  different  materials 
were  tried,  and  from  comparisons  of  these  simulations  with  corresponding  polyhedron  approximation 
simulations,  it  was  found  that  this  model  works  best  when  the  characteristic  energy  loss  per  unit  length 
dE/dx  is  simitar  for  both  materials.  Figure  6  illustrates  several  such  cases.  Unless  otherwise  noted,  all 
simulations  assume  a  2-MeV  4He+  beam,  direct  180°  backscattering,  a  1-msr  detector  solid  angle,  and 
1-pC  total  deposited  charge;  and  the  polyhedron  approximation  simulations  are  all  the  sum  of  80  separate 
flat  surface  simulations.  The  effects  of  energy  straggling  and  detector  resolution  are  not  included. 
Figure  6a  contains  simulations  of  500-A  &  on  a  Ti  sphere  and  500-A  Ti  on  a  Zr  sphere.  The  ratio  of 
dE/dx  for  2-MeV  alphas  in  these  two  materials  is  about  1.02,  and  as  a  result  it  can  be  seen  that  the  flat 
model  and  polyhedron  approximation  simulations  agree  exactly.  Figure  6b  contains  simulations  of 
2,000-A  TiOz  on  a  Ti  sphere  and  2,000-A  Ta^  on  a  Ta  sphere.  The  densities  of  the  oxides  were 
assumed  to  be  4.2  and  8.2  g/cm3,  respectively,  and  the  Ta^^a  simulation  assumes  0.2  pC  of  charge 
deposited.  These  cases  are  more  realistic,  with  dE/dx  ratios  of  about  1.13  and  0.74,  respectively,  but  only 
slight  differences  between  the  two  models  can  be  seen. 

As  stated  previously,  one  must  look  for  fairly  extreme  cases,  where  the  surface  and  bulk  values  of 
dE/dx  differ  by  about  a  factor  of  2,  in  order  to  achieve  noticeable  differences  between  the  two  methods 
of  simulation  This  is  illustrated  in  Figure  6c,  the  cases  of  300-A  of  Os  on  Ti  spheres  and  1,000- A  of  Mg 
on  Ti  spheres,  with  dE/dx  ratios  of  about  1.99  and  0.49,  respectively.  The  Mg  on  Ti  simulation  assumes 
5  pC  of  charge  deposited.  As  can  be  seen,  the  flat  model  simulations  deviate  a  small  amount  from  the 
polyhedron  approximation  simulations.  This  is  due  to  the  approximate  nature  of  the  flat  model  method. 
The  flat  model  method  overestimates  the  yield  in  the  high  energy  part  of  the  Os  peak’s  tail  and 
underestimates  it  in  the  low-eneigy  part,  as  it  always  does  when  the  value  of  dE/dx  for  the  surface  layer 
exceeds  that  of  the  bulk  material.  The  opposite  is  true  when  dE/dx  for  the  surface  is  less,  but  once  again 
the  two  values  must  differ  by  about  a  factor  of  2  before  this  is  a  problem  (e.g..  Mg  on  Ti).  Thus,  the  flat 
model  is  applicable  to  a  wide  variety  of  sample  configurations.  To  demonstrate  this,  we  determined  the 
predicted  values  of  dE/dx  for  2-MeV  4He  in  76  nongaseous  elements.  These  values  ranged  from  about 
75  keV/pm  in  Li  to  about  775  keV/pm  in  Os  and  averaged  about  389  keV/pm,  which  is  the  value  for  Ti. 
Of  these  76  elements,  65  have  dE/dx  values  that  meet  the  requirement  that  they  lie  in  the  range  between 
50  and  200%  of  the  value  for  the  "average"  element  Ti. 
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Figure  6a.  Simulations  of  180°  backscattered  2-MeV  4He  spectra  fon  500-A  Zr  on  Ti  spheres  and  500-A 
Ti  on  Zr  spheres. 
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Figure  6b.  Simulations  of  180°  backscattered  2-MeV  4He  spectra  for  2.000-A  TiOo  on  Ti  spheres  and 
2.000-A  Ta,Oc  on  Ta  spheres. 
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3.3  Ftetp-rmining  a  Radial  Profile  Directly  From  a  Flat  Profile  Fit.  For  everyday  use  as  an  analysis 
tool,  it  would  be  far  more  convenient  to  obtain  a  spectrum  from  a  compact  of  spherical  particles, 
immediately  fit  it,  assuming  a  flat  surface,  and  then  somehow  derive  information  on  the  radial  profile 
directly  from  the  fitted  equivalent  flat  profile.  In  such  a  procedure,  current  flat  surface  fitting  methods 
could  be  used  without  special  treatment  until  after  the  final  profile  is  obtained.  A  method  of  deriving  a 
radial  profile  N(r)  from  a  fitted  equivalent  flat  profile  N(0)  is  possible  for  cases  involving  spheres  with 
R  »  {.  The  following  is  a  derivation  of  such  a  method. 

When  R  »  C,  Equation  13  takes  on  the  form: 

Njfl)  =  f*  NiOOFsO-.Odr.  (25) 

'O 
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Taking  the  derivative  of  each  side  with  respect  to  6  results  in: 


3Ni 

~3f 


=  £Ni(r)^Ldr  +  [Ni(r)Fs(r,l) 


r  =  $ 


(26) 


Now,  from  Equation  15,  we  see  that  Fs  =  2 r/(P  and  thus  dFJdi  =  -4r/C3  =  -(2/{)Fs,  so  that  the  above 
can  be  rearranged  to  form: 


Nj(r)  = 


r  o  dNi 


v2y 


de 


+  Nj(4) 


t  =  r 


(27) 


Equation  27  can  be  used  to  derive  a  radial  profile  from  the  fitted  flat  profile  for  spheres  with  R  »  6, 
insofar  as  the  fitted  equivalent  flat  profile  Nj(f)  and  the  average  concentration  as  a  function  of  depth  Nj(l) 
are  tire  same.  This  requires  that  the  detector  be  at  180°  and  that  the  surface  and  bulk  values  of  dE/dx  be 
similar  (within  a  factor  of  2).  This  type  of  arrangement  is  not  possible  with  cylinders  or  when  R  -  6 
because  it  relies  on  the  fact  that  dF  Jdl  can  be  expressed  in  the  form  f({)Fs  (r,{)  (here  f (6)  =  -2/6),  so 
that  the  integral  in  Equation  26  can  be  replaced  with  f({)Nj  (C). 

The  above  method  has  practical  limitations  due  to  the  nature  of  a  typical  fitting  process.  The  first  is 
that,  using  ion  beam  techniques,  the  final  fitted  profile  is  given  in  terms  of  areal  density  A  and  atomic 
fractions  fj,  rather  than  depth  6  and  atomic  densities  Nj({).  This  problem  is  alleviated  if  we  further 
restrict  to  cases  where  surface  and  bulk  atomic  densities  are  similar  so  that  we  may  assume  that  A  a  { 
and  fj  a  Nj.  The  second  limitation  is  that  an  extremely  precise  fit  is  required  in  order  to  obtain  accurate 
information  on  the  derivative  of  a  profile  (3Nj  fd6). 

In  order  to  illustrate  the  use  of  this  method,  we  used  the  polyhedron  approximation  to  simulate  a 
spectrum  for  a  500-A-wide  linear  profile  of  Hf  on  Ge  for  the  purposes  of  fitting  with  an  equivalent  flat 
profile  and  directly  deriving  the  radial  profile  using  Equation  27.  Hf  and  Ge  were  chosen  because  they 
have  identical  atomic  densities  (N  =  4.429  x  lO^/cm2),  although  dE/dx  for  2-MeV  4He  differs  by  a  factor 
of  about  1.44  in  these  substances.  Figure  7a  shows  the  spectrum  generated  with  the  polyhedron 
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approximation  and  the  spectrum  (which  overlaps  the  first  exactly)  representing  the  fit,  assuming  a  flat 
surface.  Also  included  is  a  simulation  generated  using  the  flat  model  method  described  in  the  previous 
section,  which  deviates  only  slightly  from  the  polyhedron  simulation.  Figure  7b  shows  part  of  the  fitted 
profile  obtained  under  the  assumption  that  the  polyhedron  approximation  spectrum  was  from  a  flat  surface 
and  the  radial  profile  derived  from  it  using  Equation  27.  For  comparison,  the  flat  model  predicted 
equivalent  flat  profile  and  the  actual  radial  profile  are  also  shown.  As  can  be  seen,  the  derived  radial 
profile  occasionally  deviates  from  the  actual  profile  slightly.  Most  of  the  short-term  fluctuations  are  due 
to  the  difficulty  in  accurately  determining  the  local  values  of  the  first  derivative  of  the  profile,  which  is 
particularly  difficult  where  the  second  derivative  is  appreciable.  This  occurs  in  regions  where  there  is  a 
discontinuity  in  the  radial  profile.  Away  from  such  regions,  this  method  works  well  and,  in  fact,  the 
derived  radial  profile  effectively  agrees  with  the  actual  profile  out  to  an  areal  density  of  6  x  1018  at./cm2. 
Perhaps  somewhat  fortuitously,  the  total  amounts  of  Hf  represented  by  the  derived  and  actual  radial 
profiles  between  0  and  1  x  1018/cm2  agree  to  within  0.1%. 

In  order  to  illustrate  the  errors  that  arise  when  this  method  is  used  for  an  inappropriate  sample 
configuration,  the  above  procedure  was  performed  for  the  case  of  a  1,000-A-thick  uniform  layer  of  Ni  on 
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a  W  sphere.  While  the  atomic  density  of  Ni  (9.126  x  Kp/cm2)  and  the  atomic  density  of  W  (6.320  x 
1022/cm2)  differ  by  a  factor  of  1.44,  it  was  still  assumed  that  A  a  {  and  fj  a  Nj(fi).  The  ratio  of  dE/dx 
for  Ni  and  W  is  about  0.92,  and  therefore  this  mismatch  is  not  an  issue.  As  can  be  seen  in  Figure  8,  the 
derived  radial  profile  exhibits  larger  errors  than  before.  The  "tail"  seen  near  the  deeper  side  of  the  uniform 
layer  is  partially  due  to  the  difficulty  in  determining  dN/9 6  near  die  large  discontinuity.  However,  it  and 
the  long-term  fluctuations  seen  at  greater  depths  are  also  consequences  of  the  atomic  density  mismatch 
between  Ni  and  W.  Interestingly,  the  total  amounts  of  Ni  represented  by  the  derived  and  actual  radial 
profiles  between  0  and  1  x  1019/cm2  still  agree  to  within  0.2%. 


4.  CONCLUSION 


We  have  discussed  two  methods  of  adapting  ion  beam  techniques  for  use  on  spheres  and  cylinders, 
each  having  its  own  strengths  and  weaknesses.  Both  are  suitable  for  use  with  a  variety  of  backscattering 
techniques. 
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Figure  8a.  As  in  Figure  7,  three  simulated  Spectra  for  the  case  of  1.000-A-wide  uniform  layer  of  Ni  on 
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The  polyhedron  approximation  is  a  straightforward  and  versatile  method  as  it  puts  no  restrictions  on 
sample  composition  or  scattering  angle.  It  is,  however,  restricted  for  use  when  the  sample  radii  greatly 
exceed  the  probe  depth.  It  also  tends  to  be  quite  slow  due  to  the  large  number  of  flat  surface  spectrum 
simulations  required  in  order  to  construct  one  curved  surface  simulation.  Although  it  is  not  necessary, 
use  of  180°  scattering  proves  to  be  beneficial  as  azimuthal  symmetry  in  the  scattering  and  energy  loss 
processes  greatly  reduces  the  number  of  calculations  required.  Future  work  on  this  method  will  center 
on  speeding  the  simulation  process. 

The  first  flat  model  procedure,  by  comparison,  is  much  faster,  as  it  requires  only  one  flat  surface 
spectrum  simulation  procedure  to  be  performed  in  order  to  obtain  a  spectrum  representative  of  a  curved 
surface.  It  is  also  not  restricted  to  large  radii,  although  the  calculations  are  simpler  if  this  is  the  case.  The 
flat  model  method  is  restricted  to  a  180°  scattering  angle— a  requirement  that  can  nearly  be  met  ihrough 
the  use  of  an  annular  detector.  It  may  be  possible  to  modify  the  flat  model  so  as  to  relax  this  restriction. 
For  instance,  simply  ignoring  regions  out  of  direct  view  of  the  detector  may  allow  use  of  other  detector 
angles  in  many  cases. 

Ideally,  the  surface  and  bulk  values  of  dE/dx  should  not  differ  greatly  where  the  flat  model  is  applied. 
However,  the  level  of  performance  of  this  method  is  relatively  insensitive  to  mismatches  in  dE/dx;  the 
values  may  vary  by  as  much  as  a  factor  of  2  before  the  effect  is  noticeable.  This  is  fortunate,  in  that 
models  taking  dE/dx  variations  into  account  would  necessarily  be  much  more  complex.  When  dE/dx 
varies  between  surface  and  bulk  compositions,  then  the  locus  of  scattering  points  in  the  target  that  result 
in  backscattered  ions  with  the  same  energy  Ej  can  no  longer  be  described  by  a  simple  spherical  (or 
cylindrical)  surface  at  some  depth  C.  It  will  instead  be  a  surface  with  a  complex  shape  that  is  highly 
dependent  on  the  energy  E,  variations  in  composition,  and  the  type  of  ion  beam  used.  The  backscattering 
yield  between  energies  Ej  and  Ej  +  dE  would  then  be  proportional  to  the  number  of  atoms  contained 
between  two  such  complex  surfaces.  Thus,  it  would  be  more  appropriate  to  assign  to  each  sublayer  in  an 
equivalent  flat  profile  the  same  average  composition  as  found  between  two  such  "surfaces  of  constant  ion 
energy,"  rather  than  between  two  "surfaces  of  constant  depth"  as  described  here.  The  observed 
insensitivity  of  this  purely  geometric  method  to  dE/dx  variations  implies  that  exact  descriptions  of  the 
shapes  of  these  surfaces  are  not  crucial.  Calculating  these  shapes  and  the  resultant  yields  would  be  an 
enormous  task  and  far  more  difficult  than  a  flat  surface  spectrum  simulation.  It  would  also  be  difficult 
to  set  up  a  computer  program  to  do  this  for  multiple  types  of  ion  beam  analysis  and  sample  configurations. 
In  contrast,  the  flat  model  method  described  here  amounts  to  the  addition  of  a  single,  simple  step  (the 
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construction  of  an  equivalent  flat  surface  profile)  to  current  fiat  surface  spectra  fitting  procedures  and 
generally  does  not  even  require  reprogramming. 

The  second  flat  model  method,  for  deriving  a  radial  profile  directly  from  a  fit  obtained  under  the 
(incorrect)  assumption  that  the  surface  is  flat,  is  at  present  further  restricted  to  situations  where  the  atomic 
density  does  not  vary  greatly  in  the  sample  and  is  also  restricted  to  large  radii  spheres.  However,  it  may 
be  the  fastest  and  most  convenient  method  discussed  here,  and  it  would  yield  useful  qualitative  information 
even  in  situations  where  it  is  not  entirely  appropriate.  It  would  be  particularly  well  suited  to  resonant 
methods  such  as  the  hydrogen  profiling  technique  based  on  the  H(15N,ay)12  C  nuclear  reaction  (Horn  and 
Lanford  1988),  where  the  observed  profile  is  obtained  directly;  no  spectra  need  be  simulated.  Also,  if  the 
detected  particle  is  a  y  ray,  there  would  be  no  restriction  on  scattering  angle.  In  the  future,  we  hope  to 
adapt  a  number  of  different  ion  beam  techniques  for  use  in  the  study  of  powder  surface  contamination. 

As  not  all  powder  particles  are  spherical  or  cylindrical,  it  would  also  prove  useful  to  consider  the  case 
of  randomly  shaped  powder  particles.  Our  efforts  to  date  suggest  that  it  may  not  be  necessary  to  know 
the  actual  shapes  of  the  particles,  so  long  as  the  distribution  of  surface  normal  angles  presented  to  the 
beam  is  known.  It  may  well  be  that  a  surface  consisting  of  a  large  number  of  randomly  shaped  and 
randomly  oriented  powder  particles  will  tend  to  have  the  same  distribution  of  normal  angles  as  a  sphere, 
so  that  spherical  surface  analysis  techniques  can  be  used  without  modification.  If  large  and  small  surface- 
normal-to-beam  angles  are  more  equally  represented,  then  cylindrical  surface  analysis  techniques  might 
be  appropriate.  We  intend  to  continue  to  explore  the  relationship  between  particle  shape  and  surface 
normal  angle  distributions,  so  that  spectra  from  other  powder  particle  shapes  may  be  simulated. 
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